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ABSTRACT 

The 10/im silicate feature observed with Spitzer in active galactic nuclei (AGN) reveals some puzzling be- 
havior. It (1) has been detected in emission in type 2 sources, (2) shows broad, flat-topped emission peaks 
shifted toward long wavelengths in several type 1 sources, and (3) is not seen in deep absorption in any source 
observed so far. We solve all three puzzles with our clumpy dust radiative transfer formalism. Addressing (1), 
we present the spectral energy distribution (SED) of SST1721H-6012, the first type 2 quasar observed to show 
a clear lO/im silicate feature in emission. Such emission arises in models of the AGN torus easily when its 
clumpy nature is taken into account. We constructed a large database of clumpy torus models and performed 
extensive fitting of the observed SED. We find that the cloud radial distribution varies as r"^^ and the torus con- 
tains 2-4 clouds along radial equatorial rays, each with optical depth at visual ^60-80. The source bolometric 
luminosity is ^3 • lO'^L©. Our modeling suggests that < 35% of objects with tori sharing these characteris- 
tics and geometry would have their central engines obscured. This relatively low obscuration probability can 
explain the clear appearance of the 10/im emission feature in SST1721H-6012 together with its rarity among 
other QS02. Investigating (2) we also fitted the SED of PG1211H-143, one of the first type 1 QSOs with a 
10/im silicate feature detected in emission. Together with other similar sources, this QSO appears to display an 
unusually broadened feature whose peak is shifted toward longer wavelengths. Although this led to suggestions 
of non-standard dust chemistry in these sources, our analysis fits such SEDs with standard galactic dust; the 
apparent peak shifts arise from simple radiative transfer effects. Regarding (3) we find additionally that the dis- 
tribution of silicate feature strengths among clumpy torus models closely resembles the observed distribution, 
and the feature never occurs deeply absorbed. Comparing such distributions in several AGN samples we also 
show that the silicate emission feature becomes stronger in the transition from Seyfert to quasar luminosities. 
Subject headings: dust, extinction — infrared: general — galaxies: active — quasars: individual: PG1211H-143 
— quasars: individual: SST1721H-6012 — radiative transfer — 



1. INTRODUCTION 

Unified schemes of active galactic nuclei (AGN) require an 
obscuring dusty torus around the central source, giving rise to 
a type 1 line spectrum when there is direct view of the cen- 
tral engine and type 2 characteristics when it is blocked (e.g. 
Antonucci 1993; Urry & Padovani 1995). The torus, which 
is comprised of dusty clouds that are individually optically 
thick (Krolik & Begelman 1988), reprocesses the radiation 
it absorbs into longer wavelengths, creating a distinct signa- 
ture in the observed infrared. Silicates, a major constituent 
of astronomical dust, reveal their presence through the spec- 
tral feature at 10 /im. Among type 1 AGN, QSOs display 
the feature in emission (Siebenmorgen et al. 2005; Hao et al. 
2005; Sturm et al. 2005), while average SEDs of Seyfert 1 
galaxies have either a flat 10/im feature (Wu et al. 2009) or 
show it in mild absorption (Hao et al. 2007). Seyfert 2 galax- 
ies generally display an absorption feature with limited depth, 
much shallower than in ultra-luminous IR galaxies (e.g. Hao 
et al. 2007; Levenson et al. 2007). An intriguing result comes 
from the Spitzer observations of seven high-luminosity type 2 
QSOs by Sturm et al. (2006). While individual spectra appear 
featureless, the sample average spectrum shows the 10/im fea- 
ture in emission. 

Heated dust will produce the feature in emission whenever 
it is optically thin. When the dust optical depth at 10 /im ex- 
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ceeds unity, the feature still appears in emission in viewing 
of the illuminated face of the dust but in absorption when the 
dust is between the observer and heating source. In the ab- 
sence of a formalism for radiative transfer in clumpy media, 
early models of the AGN torus employed smooth density dis- 
tributions instead (e.g. Pier & KroHk 1992, 1993; Granato & 
Danese 1994; Efstathiou & Rowan-Robinson 1995; Granato 
et al. 1997; Fritz et al. 2006). These models predict that type 1 
sources, where the observer has a direct view of the torus in- 
ner, heated face, will generally produce an emission feature, 
although some examples of absorption features do exist (Pier 
& Krolik 1992; Efstathiou & Rowan-Robinson 1995). Type 2 
viewing in most cases produces an absorption feature, whose 
depth is quite large on occasion, much larger than ever ob- 
served. An emission feature is rarely produced from such 
viewing, but it should be noted that it is under certain condi- 
tions for instance in the model presented by Fritz et al. (2006). 
A formalism for handling clumpy media was developed by 
Nenkova et al. (2002, 2008a) (hereafter N02 & N08a); the 
formalism holds for volume filling factors as large as 10%. 
Their models show that a clumpy torus will never produce a 
very deep absorption feature and that the feature displays a 
much richer behavior than in smooth density models; in par- 
ticular, type 1 viewing can produce an absorption feature in 
certain models and type 2 viewing can lead to an emission 
feature in others (Nenkova et al. 2008b, N08b henceforth). 

While the Sturm et al. (2006) data suggest the possibility 
of a 10/im emission feature in QS02, the only unambiguous 
evidence for such a feature in a type 2 AGN was presented re- 
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cently for the Seyfert galaxy NGC 21 10 (Mason et al. 2009).^ 
Here we present the first unambiguous case of an emission 
feature in a type 2 quasar, SST1721+6012, and perform ex- 
tensive fitting of its spectral energy distribution (SED) with 
clumpy torus models. 

The comparison of torus model predictions with observa- 
tions is somewhat problematic because the overwhelming ma- 
jority of these observations do not properly isolate the torus 
IR emission. Starburst emission is a well known contaminant 
in many cases, and we selected SST1721H-6012 for modeling 
precisely for this reason as its spectrum seems free of starburst 
indicators. However, even IR from the immediate vicinity of 
the AGN may not always originate exclusively from the torus. 
High-resolution observations of NGC1068 by Cameron et al. 
(1993) and recently by Mason et al. (2006) demonstrate that 
the torus contributes less than 30% of the IQfim flux collected 
with apertures > 1" in this object, with the bulk of this flux 
coming from dust in the ionization cones (Braatz et al. (1993) 
also found that at least 40% of the 12.4/im flux in this source 
do not originate from the torus). The significance of IR emis- 
sion from the narrow line region (NLR) was noted also by 
Schweitzer et al. (2008). However, because the dust in the 
ionization cones is optically thin, its IR emission is isotropic 
and does not generate differences between types 1 and 2. Ob- 
servations show that such differences do exist. In particular, 
the Hao et al. (2007) compilation of Spitzer IR observations 
shows a markedly different behavior for the IQ^m feature be- 
tween Seyfert 1 and 2 galaxies. Accepting the framework 
of the unification scheme, these differences can be attributed 
only to the torus contribution. Thus it seems that, unfortu- 
nately, a general rule does not exist and the situation must be 
investigated case by case. Our aim here is to examine whether 
the torus contribution alone can reproduce the observed SED 
of SST1721H-6012, yielding a range of possible parameter val- 
ues that describe the dusty cloud distribution in this source 
(§2). 

We also investigate the cause for apparent shifts of the sili- 
cate feature peaks towards long wavelengths (§3). Such shifts 
have been reported for sources that show the lO/im feature 
in emission (Siebenmorgen et al. 2005; Sturm et al. 2005; 
Hao et al. 2005), and attributed to non-standard dust chem- 
istry. However, these shifts were never seen in absorption, 
suggestive of radiative transfer effects instead. Finally, in 
§4 we compare the observed distribution of silicate feature 
strengths among the Hao et al. (2007) sample of AGN with 
the synthetic distribution of feature strengths in our database 
of clumpy torus model SEDs. 

2. SILICATE 10-MICRON EMISSION FEATURE IN QS02 

Although not expected in type 2 sources, possible detection 
of the lOfim emission feature was reported by Sturm et al. 
(2006). The feature was only identified after averaging the 
SEDs of a number of type 2 QSOs, which individually show 
no significant indication of the feature. Recently Mason et al. 
(2009) presented the first unequivocal detection of an emis- 
sion feature in an individual type 2 source, the Seyfert galaxy 
NGC 2110. We present the Spitzer SED of the type 2 quasar 
SST1721H-6012 that shows the 10/im and 18/im silicate fea- 

^ Teplitz et al. (2006) have suggested a lOfim emission feature in the 
Spitzer spectrum of QS02 FSC10214+4724. The suggestion is problematic 
because the object's redshift is so high (z = 2.2856) that the iOfim feature was 
not fully in the spectral range of the IRS instrument. The rest-frame spectrum 
is cut off around 12/^m, before the continuum longward of the feature could 
be established. 



tures in emission. In this section we report on the results of 
fitting the SED of SST1721H-6012 with clumpy torus models, 
and derive multiple parameters characterizing the source. 

2.1. Observations 

The source SSTXFLS J172123.1H-601214 was first identi- 
fied as an AGN candidate in the Spitzer First Look Survey 
(FLS) by Lacy et al. (2004). It has a redshift of z = 0.325, 
and was not present in the SDSS at that time. In 2007, Lacy 
et al. (2007a) categorized it as a type 2 quasar based on the 
presence of optical, narrow [N V] emission lines and through 
emission line ratio diagnostics introduced by Baldwin et al. 
(1981). In the same year Lacy et al. (2007b) presented, to- 
gether with other sources, a wide-range SED for this source, 
including a mid-IR spectrum taken by the Infrared Spectro- 
graph (IRS) aboard Spitzer. 

The IRS observations (Astronomical Observation Request 
1406768) were taken on 2005 August 14 in staring mode us- 
ing the short and long low resolution modules to obtain con- 
tinuous coverage from 5.2-38 /im, and were passed through 
the S14.0 version of the SSC pipeline. The signal-to-noise 
ratio varied through the spectrum, the deepest observations 
being targeted on the redshifted wavelengths of the strong 
spectral features expected to lie in the 7-15 /im range. In 
the short wavelength module, two 14s ramps were taken in 
second order, and a single 60s ramp in first order In the long 
wavelength module, two 30s ramps were taken in both first 
and second order. The spectra from each module were op- 
timally extracted using SPICE.'* The resulting spectra were 
trimmed, combined, and resampled in constant energy bins of 
AA/A ~ 0.01, resulting in a spectrum ranging from 4.0 
to 27.1 /im (rest wavelength). Uncertainty estimates from 
SPICE were propagated through the process in the usual man- 
ner. For the fitting we excluded a few data points at shorter 
and at longer wavelengths due to poor signal-to-noise ratio. 
Additionally, we make use of two photometric data points 
from the Infrared Array Camera (IRAC) component of the 
Spitzer First Look Survey (Lacy et al. 2005) at rest wave- 
lengths of 2.7 /im and 3.4 /im, both with very small intrinsic 
uncertainties, as they greatly help defining the shape of the 
SED in the regions of hot dust emission. Cross-calibration be- 
tween IRS and IRAC is accurate to better than 10% (L. Yan, 
personal communication). 

Despite a certain noisiness in the IRS spectrum, a clear 
presence of polycyclic aromatic hydrocarbon features (PAHs) 
can be safely excluded. Considering additionally its lack of a 
[Ne 11] emission line at 12.8 /im, the spectrum shows no signs 
of star formation. Furthermore, the SED seems free of other 
emission lines, with one disputable exception. Locally, the 
flux peaks around 10.5 /im, which coincides with the [S IV] 
emission line at 10.51 /im reported to be found in 11 out of 
12 type 2 sources by Zakamska et al. (2008). This radiation, 
if indeed credited with an emission line, would stem from the 
AGN itself, but our spectrum does not show any other lines 
originating from the AGN, like [Ne III] at 15.5 /im and [Ne V] 
at 14.3 /im. Within the frame of this work, we therefore at- 
tribute the peak flux at ^10.5 /tm entirely to silicate emission. 

2.2. Modeling 

N02 & N08a describe an analytic formulation of radia- 
tive transfer in a clumpy, dusty medium heated by a radi- 

^ http://ssc.spitzer.caltech.edu/postbcd/spice. 
html 
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TABLE 1 

Clumpy parameters used in fitting 



Parameter 


Sampled Values 




0, 0.5, 1, 1.5,2,2.5,3 


No 


1-25 


TV 


5, 10, 20, 30, 40, 60, 80, 100, 150, 200, 300, 500 


a 


15, 20, 25, 30, 35, 40, 45, 50, 55, 60, 65, 70, 75, 80 


i 


0, 10, 20, 30, 40, 50, 60, 70, 80, 90 


Y 


2 - 5, 10, 20, 30, 40, 50, 60, 70, 80, 90, 100, 150, 200 



Note. — a and i are measured in degrees. 

ation source. The formalism was implemented in the code 
Clumpy, which takes as input a toroidal distribution of point- 
like dust clouds around a central source. The dust in each in- 
dividual cloud has an optical depth Ty, defined at 0.55 /im, 
and standard ISM composition of 47% graphite with optical 
constants from Draine (2003) and 53% "cold" silicates from 
Ossenkopf, Henning, & Mathis (1992) (OHM hereafter). The 
dust sublimation temperature defines the torus inner radius Rd 
and is set to 1500 K. The cloud distribution is parametrized 
with the radial power law l/r^ between R^i and the outer ra- 
dius YRd, where q and Y are free parameters. Another free 
parameter is A^o, the average number of clouds along a radial 
equatorial ray. In polar direction the number of clouds per 
radial ray is characterized by a Gaussian, so that at angle (3 
from the equatorial plane it is A^o e"(''/°') , with a the last free 
parameter of the cloud distribution. The final parameter is /, 
the observer's viewing angle measured from the torus axis. 

We employed Clumpy to produce a large database^ of 
model SEDs f\ = \F\/Facn, with Facn the total bolometric 
flux. The observations provide a set of fluxes, F", at wave- 
lengths Xj, j = I ■■ -N. Our fitting procedure involves search- 
ing the entire database for the model that minimizes the error 
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where Aj are individual errors on the XjF" 



and f" are the 

model fluxes at the same set of wavelengths as the data. Each 
model SED is scaled by the factor /^icw that minimizes E, 
determining the AGN bolometric flux for this model. Since 
the data dynamic range is only w 3, the fitting procedure can 
be safely executed in linear space. 

2.3. Results 

We calculated E for all the CLUMPY models whose param- 
eters are listed in Table 1, resulting in a database of more than 
4.7 million entries. This large set contains as a subset all the 
parameters that N08b found to be plausible. Figure 1 shows 
the data and the best-fitting CLUMPY model. The two photo- 
metric IRAC points play a crucial role in the fits by expanding 
the data into the short wavelengths. 

Although the model presented in Figure 1 produces the 
smallest nominal error E, a number of other models have er- 
rors that differ from it only in the third significant digit. Be- 
cause of the large degeneracy of the radiative transfer problem 
for heated dust, the SED is a poor constraint on the properties 
of the source; a meaningful determination of model param- 
eters requires also high-resolution imaging at various wave- 
lengths (e.g., Vinkovic et al. 2003). The axially symmetric 

' Models are available at http ; / /www . pa . uky . edu/ clumpy 
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Fig. 1 .— SED of SST1721+6012. Spitzer IRS data are shown in dark gray 
with the errors in light gray shade. Two IRAC photometry points are marked 
with crosses. The black line shows the best-fit CLUMPY model, which pro- 
duces an error £„,|„ = 0.212 (see Equation 1). Its parameters are q = 1.5, 
Nq = 3, TV = 80, cr = 20, y = 30 and i = 60. The inset shows the data and the 
best fit model using \Fx and linear scales for a better display of the iOfim 
emission feature. 



clumpy torus model requires a relatively large number of in- 
put parameters, further exacerbating the degeneracy problem. 
We define E,- = 100 • (E - E^nm) / E„iin as the relative deviation 
of a model from the best-fit one. Then, 199 models have 
Ef < 5%, within a fraction of the minimal error E„^i„ = 0.212, 
and the bar diagrams of these models are shown in the top 
rows of Figures 2 and 3 for each of the six parameters. All 
but two of these models share the same value of q = 1 .5, indi- 
cating that this parameter can be considered well constrained. 
Similarly, for 90% of all models No is either 3 or 4, so this 
parameter is only slightly less well constrained. The distribu- 
tions of the parameters Ty, a and / are broader, but still show 
well defined peaks. For these parameters we can only deduce 
a plausible range. In contrast, the parameter Y has a flat distri- 
bution that covers every sampled value Y > 20; this parameter 
is undetermined, except for the indication of a lower bound. 

The choice E,- < 5% is of course arbitrary. Increasing 
slightly the range of accepted models, the bar diagrams can 
be expected to remain peaked if the parameters are well con- 
strained. The figures show that this is indeed the case for q and 
A'o, whose distributions remain reasonably peaked in spite of 
the large increase in the number of accepted models (almost 
13,000 are selected by the criterion Er < 20%). To a lesser 
degree, this is also the case for Ty. In contrast, the distribu- 
tions of a and i, which also start out peaked, flatten out signif- 
icantly as the acceptance criterion is relaxed, indicating that 
SED analysis lacks the predictive power to constrain these 
parameters in SST1721-I-6012. The only meaningful results 
are that, in all likelihood, a < 50° and / < 70°, i.e., edge-on 
viewing is excluded. Furthermore, these parameters are not 
entirely independent of each other since a clear line of sight 
to the AGN can be obtained for different combinations of the 
two. In fact, the interdependence of a, i, and to some de- 
gree A^o constitutes the greatest source of degeneracy within 
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TABLE 2 

Properties of fitted parameters for SST1721+6012 



1 2 
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Fig. 2. — Bar diagrams of thi'ee CLUMPY parameters well constrained by 
fitting. From left to right, the columns correspond to q. No, and tv. The 
parameters were sampled as listed in Table 1 . Rows correspond, from top to 
bottom, to an increasing acceptance on the fitting en'or relative to the best-fit 
model, as marked on the right, with the resulting number of models increasing 
accordingly — 199, 1691, 5210 and 12854. The height of the bar at any value 
of a parameter is the fraction of all accepted models. 



1.0 



0.5 



1.0 



1.0 



0.5 



1.0 



0.5 



0.0 



-III ■ 


,,,iilL,, 


llllllllffllllflllll 


III 


M 1 1 1 1 1 ^ 


lIlMIIMIiiiiliiiil 


II 


1 m 1 1 1 1 , , 


MMMIMIiiiiIiiiiI 


Ill 


M 1 1 1 M 1 1 , 


llMIIIMIiiiiliiiil 



15 30 45 60 75 30 60 90 50 100 150 200 

o- [deg] i [deg] Y 

Fig. 3 . — Same as Figure 2, but for the three less well constrained CLUMPY 
parameters cr, i, and Y. These distributions flatten out more quickly with 
growing acceptance error. 

the clumpy torus SEDs. The final parameter, the torus radial 
thickness Y, is undetermined. As noted already in N08b, the 
SEDs of models with a steep radial cloud distribution > 1) 
are insensitive to increasing Y because most of the clouds are 
concentrated in the torus inner region. The only constraint we 
can deduce is the lower bound F > 10, indicating that the torus 
could be compact, in agreement with other AGN observations 
(see N08b and references therein). 

Table 2 summarizes the likely values constrained by fitting. 
We cannot give exact confidence intervals since our distri- 
butions are not continuous. If a parameter is perfectly con- 
strained, all models then have the same value. Denoting by 
H the fraction of models at the distribution peak value, such 
a parameter would have H = 1. On the other hand, a flat dis- 
tribution over the entire range of sampled values indicates a 
completely non-constrained parameter. If the number of pa- 
rameter values in the sampled range is B, the height of each 
bar would then be 1 /B. Introduce w = H/B. A perfectly con- 
strained parameter will have w = I {= H = B) while for an 



Parameter 


Best Fit 


Peak" 


90%-Range'' 


B'^ 


H<1 


w^ 


q 


1.5 


1.5 


1 - 1.5 


2 


0.81 


0.40 


No 


3 


3 


2-4 


3 


0.42 


0.14 


TV 


80 


80 


30-100 


5 


0.30 


0.06 


¥ 


30 


20 


20 - 200 


11 


0.13 


0.01 




20 


15 


15-40 


6 


0.27 


0.05 


i 


60 


50 


0-60 


7 


0.16 


0.02 



Note. — Statistical indicators for the sample of all 1691 models with 
E,- < 10% deviation from the best-fit model. 
" Value of the parameter at the distribution peak. 

Range around the peak containing at least 90% of the sample models. 

Number of sampled values in the 90%-range. 

Fraction of all accepted models at the distribution peak. 

Measure of how well the parameter is constrained (see text); the closer w is 
to unity the higher is the significance of the determined value. The values for 
the last two entries cannot be directly compared to all others since the range 
of both cr and (' is finite whereas for all other pai'ameters it is in principle 
unhmited. 

unconstrained parameter w = 1 /B^, decreasing when the num- 
ber of sampled values is increasing. We select as our sample 
the 1691 models with Er < 10%. While admittedly arbitrary, 
this selection ensures a strict acceptance criterion while still 
giving a statistically large sample. For each of the model pa- 
rameters we identify the minimal interval around the distri- 
bution peak containing at least 90% of the sample's models. 
These ranges are listed in Table 2, together with the number 
of sampled values (bars) in these intervals, which is our mea- 
sure of B. The last column lists the corresponding values of 
w, reinforcing the perception conveyed by the bar diagrams 
regarding the degree of confidence (or lack thereof) in each of 
the derived model parameters. 

This analysis shows that the radial cloud distribution in 
SST1721H-6012 is well constrained at <7 = 1.5; although the 
90% range contains also q = 1, 81% of the models are at 
q= 1 .5. The likely value of A^o is similarly well constrained to 
the range 2-4; even though this parameter was densely sam- 
pled in steps of 1 all the way to 25, half of the 50,000 best 
models fall within this narrow range. The third reasonably 
well determined parameter is Ty ~ 80, whose likely value is 
between 30 and 100. Note that the values of these parameters 
for the best-fitting model are <7 = 1.5, A^o = 3 and Ty = 80, and 
that the close agreement with the distribution peaks is not a 
given — in principle, the best-fit model could fall anywhere 
inside the acceptable ranges. We have tried to put stronger 
constraints on the less well-defined parameters a, i, and Y, 
by holding the values of the relatively well-constrained pa- 
rameters fixed at ^ = 1.5, No = 2-4, and Ty = 60- 100. This 
had little effect on the distributions of the unconstrained pa- 
rameters, although the a bar-diagrams became slightly more 
peaked, showing a hint of greater preference for cr w 15-30. 
We conclude that it is impossible to deduce cr, /, and Y for 
SST1721H-6012 from SED analysis alone. 

2.4. Source Type 

In the standard form of the unification approach, the clas- 
sification of an AGN as type 1 or 2 is uniquely determined 
by the relation between the viewing angle / and the torus an- 
gular thickness a. In a clumpy medium, on the other hand, 
the source type is a matter of probability. Denote by A^(/) the 
average number of clouds along a radial ray at angle /, then 



^esc(0 = e 



(2) 
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Fig. 4. — Histograms of P^hx^ the probability that the AGN is obscured 
by the torus in the CLUMPY model (the probability that the model produces 
a type 2 source). Each panel corresponds to a different maximal acceptance 
error, as marked at the top together with the corresponding number of models. 
Each bin width is 0.1 and its height is the fraction of all accepted models. The 
mean value in each panel, from left to right, is 0.10, 0.1 1, 0.17 and 0.22 

is the probability that a photon emitted by the AGN will es- 
cape the torus. The source has a probability Pe,sc(0 to appear 
as a type 1 AGN and Pobsc(0 = 1 -^esc(0 as a type 2. With our 
Gaussian parametrization for the cloud angular distribution, 

A^(0 = M,e-t(^°-''/-l\ (3) 

The AGN type is probabilistic, and it depends on /, a and A^o- 
Since SST1721+6012 is a type 2 quasar, the a priori ex- 
pectation would be that Pobsc is large. We find this not to be 
the case. The best-fit model has Pobsc = 27%, and more than 
75% of all models with E, < 15% have Pobsc < 33%. Fig- 
ure 4 displays the histograms of Pob,sc for the models accepted 
at various tolerance levels, showing that the majority of mod- 
els have fobsc < 10% (in the first 3 panels). Such low proba- 
bility would pose a problem if these were the numbers for a 
large sample of type 2 sources. However, SST1721-I-6012 is a 
relatively rare type 2 quasar with a clear lOfim emission fea- 
ture; of the more than twenty QS02 with measured IR SEDs, 
NGC 2110 is the only other source with such unambiguous 
emission feature. The emission feature requires a direct line 
of sight to a significant fraction of the hot surfaces of directly 
illuminated clouds on the far inner side of the torus. Because 
obscuration of the AGN involves a single line of sight while 
the IR flux measurements integrate over many lines of sight, 
the relatively low values of Pobsc that emerge from the model- 
ing are commensurate with the clear appearance of the 10/^m 
emission feature in SST1721-I-6012 and its rarity among other 
QS02 (see also §3). 

2.5. AGN Luminosity 

Since the central engine is obscured in SST1721-I-6012, a 
direct measurement of the AGN bolometric luminosity is im- 
possible. However, the bolometric flux enters directly into 
the fitting procedure (see Equation 1) as the scale factor that 
minimizes the error in matching the model spectral shape 
with the data. The source luminosity Lacw is then derived 
from its luminosity distance Dl = 1.703 Gpc, obtained from 
the redshift z = 0.325 for standard cosmological parameters 
(Ho = 70 km s"' Mpc"', = 0.3, flat universe). The best-fit 
model has L12 = Lagw/IO'^L© = 3.47, and Figure 5 shows the 
distribution of logL^ derived for all fitted models within a 
given acceptance error. At the most restrictive level, all the 
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Fig. 5. — Histograms of the logarithm of the AGN bolometric luminosity. 
^12 = ^agn/IO'^^q. derived from the scaling of each CLUMPY model (see 
Equation 1). Each panel con'esponds to a different maximal acceptance en'or, 
as marked at the top together with the corresponding number of models. Each 
bin is 0.2 wide, and its height is the fraction of all accepted models. The mean 
value of L12 in each panel, from left to right, is 3.45, 2.71, 2.37 and 2.14. 

models fall in the range 1 . 1 < L12 < 6.5 and the mean value is 
3.45, similar to the best-fit model. As the acceptance becomes 
less restrictive, the range of accepted models extends to lu- 
minosities lower than lO'^L©, but its upper boundary stays 
unchanged; the figure panels for < 10% and E^ < 20% 
for instance are very similar except for the presence of more 
L < IO'^Lq models in the latter The reason is simple. The 
luminosity scale factor is J F\dX, and as is evident from Fig- 
ure 1, a large fraction of the integral is contained at wave- 
lengths that are missing from the data as they are shorter than 
the IRAC measurements. Model SEDs that drop precipitously 
before the IRAC points can still produce a small error esti- 
mate E by reasonably fitting all other, longer wavelengths. 
Such models will be formally acceptable-but only because the 
short wavelength region, crucial for the luminosity determina- 
tion, is so poorly sampled in the data. Observations at these 
short wavelengths will constrain better the SED, and provide 
a more accurate determination of L12. With the current data, 
our best estimate is L12 ~ 3 with a likely range of 1-7. 

Integrating the mid-infrared (MIR) luminosity only. Lacy 
et al. (2007b) find Lmir = 0.25 • lO'^L©. Richards et al. (2006) 
show that the bolometric correction from the mid-infrared is 
about a factor of eight. With this correction, the earlier es- 
timate gives L12 ^ 4, in good agreement with the detailed 
Clumpy calculations. 

3. FEATURE SHAPE AND ORIGIN OF 10-MICRON EMISSION 

After many years in which it remained undetected in type 1 
AGN, the 10/xm feature was finally discovered in emission 
in Spitzer observations (Siebenmorgen et al. 2005; Hao et al. 
2005; Sturm et al. 2005). In addition, the 18/im feature ap- 
pears in quite prominent emission. All three teams noted the 
large differences with Galactic sources — the lOfim emis- 
sion feature in AGN is much broader, and in most cases 
its peak seems to be shifted to longer wavelengths, up to 
~11 /im. Analyzing the feature with the simple approxima- 
tion K\B\{T) (optically thin emission from dust at the sin- 
gle temperature T), all three teams found significant differ- 
ences between the dust absorption coefficient in AGN and 
the interstellar medium, suggestive of a different mix of the 
silicate components. Significantly, though, the shifts toward 
longer wavelengths apparent in emission features were never 
reported in absorption; AGN absorption features reach their 
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Fig. 6. — Apparent shift of feature peak in quasar PG1211+143 as a ra- 
diative transfer effect. Top: Tlie Spitzer data (black) sliow the Wpm and 
silicate features in emission. The SED of the best-fit CLUMPY model 
(within 8-30 p,m) is shown in green; its parameters are q = {).NQ = 5,Ty = 20, 
(T = 25, Y = 20, i = 60. The model reproduces observations that prompted 
suggestions for non-standard dust composition. Two underlying continua are 
constructed as splines with (blue) and without (red) mid-range pivots over 
the 14—14.5 /jm inter-feature region. Middle: Continuum-subtracted fluxes 
for each of the continua in the top panel. Bottom: The flux ratio F jFcom for 
each continuum. 

deepest level at the same wavelengths as Galactic sources, 
^9.8 /im. The different behavior of emission and absorp- 
tion features suggests that the apparent peculiarities of AGN 
emission features do not arise from the dust composition, but 
rather from radiative transfer effects. 

To further investigate this, we analyzed the Spitzer data of 
quasar PG121 1-1-143, one of the sources in the original dis- 
covery paper of Hao et al. (2005), which shows both silicate 
features in emission. While the data cover 5-35 /im, for the 
model fitting we employed only wavelengths between 8 and 
30 [im. Fitting the shorter wavelengths with our torus mod- 
els proved rather difficult in this source. The same problem 
arises in other PG quasars, where Mor et al. (2009) find that 
high flux levels at short wavelengths necessitate the addition 
of a hot dust component to the torus emission in their models 
(see also Netzer et al. 2007). Wavelengths longer than 30 /im 
were omitted in the fitting because of high noise levels. The 
top panel of Figure 6 plots the observed SED of PG1211-I-143 
between 5 and 30 [im in black color and the SED of the best- 
fit model in green. In addition to a prominent 18/im feature, 
the displayed model shows a broad 10/im feature that reaches 
local peak emission at 11.6 /im. An analysis of the feature 
shape requires the construction of an underlying continuum. 
Sirocky et al. (2008) discuss this problem in detail and show 
that the proper continuum definition requires a spline fitted 
to the two wavelength regions shorter than the 10/im feature 
and longer than the 18/im feature, and in between them. This 
spline fit to the model results is shown with blue color in the 
figure. Although the central region, 14—14.5 /im, is essen- 
tial for a correct definition of the continuum, it was missing 



TABLE 3 

Properties of fitted parameters for PG121 1-h143 



Parameter 


Best Fit 


Peak" 


90%-Range'' 


B'^ 






y 


20 


20 


20 


1 


1.00 


1.00 


1 








0-0.5 


2 


0.61 


0.30 


TV 


20 


20 


20-30 


2 


0.50 


0.25 


No 


5 


6,7* 


2-9 


8 


0.21 


0.03 


a 


25 


25 


15-60 


10 


0.29 


0.03 


i 


60 


0,10,20,40* 


0-70 


8 


0.14 


0.02 



Note. — Statistical indicators for the sample of all 28 models with relative 
deviation E,- < 10% from the best-fit model, and with Nq < 10 (see text), 
listed in descending order of constraint. Footnotes a-f identical to Table 2. 

Peak comprises multiple bins; all listed bins have equal heights. 



from earlier analyses. The corresponding spline is plotted 
in red for comparison. The figure middle panel shows the 
continuum-subtracted flux in each case. The feature peaks 
at 10.0 /im in the properly constructed continuum, but has 
a flat plateau between ^9.8-11.6 /im that peaks nominally 
at 10.5 fim in the traditional continuum. The bottom panel 
shows the ratio F /Fcont for each continuum. Under the com- 
mon parametrization with kxBx{T), the F//^cont curves would 
be taken as the actual dust absorption coefficient. However, 
they are the outcome of radiative transfer calculations with 
the standard OHM dust, whose absorption profile looks quite 
different; these artificial "absorption coefficients" are much 
flatter than the peaked shape of the input k\. 

The reason for the peculiar shape of the emission feature 
is quite simple. The feature originates from the optically thin 
emitting layer on the bright surfaces of clouds illuminated di- 
rectly by the AGN. Absorption by other clouds encountered 
on the way out toward the observer alters the feature's shape. 
This absorption is strongest at the feature peak, where the ab- 
sorption coefficient is largest, and Ty ^ 20 is where single 
clouds become optically thick at that peak. When the gen- 
erated photons encounter ^1 cloud along the remaining part 
of the path toward the observer, the peak is absorbed while 
photons in the wings escape freely, effectively flattening the 
shape of the feature. An increasing number of clouds along 
the path would absorb the peak and the wings of the feature 
more strongly, producing a self-absorption dip in the feature's 
shape, and eventually suppressing the entire feature (see also 
Fig. 2 in N08b). The apparent shift toward longer wavelengths 
arises from the interplay with the shape of the rising contin- 
uum underneath the feature. It may be noted that such ap- 
parent variations in the shape of the silicate emission feature 
in evolved stars prompted the suggestion of dust chemical 
evolution (Little-Marenin & Litde 1990; Stencel et al. 1990), 
but were similarly shown to reflect radiative transfer effects 
(Ivezic & Elitzur 1995). 

Table 3 summarizes the analysis of the distribution of all 
Clumpy models with fitting errors within 10% of the best- 
fit model. This prescription is identical to the one employed 
for SST1721-I-6012, but in the present case it yields only 38 
models instead of 1691. These models further break into two 
distinct groups with different ranges of A^o, the radial number 
of clouds in the equatorial plane. While 28 models have A^o < 
9, the other 10 fall in the No = 16-18 range, with a large gap 
between the two groups. Because values of A'o larger than 
~10 are unlikely in general (see §3.4 in N08b), we exclude 
the ten models with > 16 from our sample. 

As before, only three parameters are well constrained. The 
radial cloud distribution is again well-constrained, but now 
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Fig. 7. — Cloud column (left) and lO^m emission (right) along viewing rays through the best-fit torus model for PG1211+143 (see Table 3). (al): Map of 
the cloud column. Axes are linear displacements xjRi and yjRd from the central AGN, with the dust sublimation radius. The gray scale is linear, with 
white standing for zero clouds and darker shades indicating higher cloud columns. {a2): One-dimensional cut through the cloud number distribution in (al) 
along the A-coordinate at y = 0. {a3): Same as (a2), but vertically along the y-coordinate at x = 0. (bl): The distribution of Wfim emission emerging from the 
central lORj X lORj. The gray scale is linear, darker shades indicating higher 10/^m emission. (b2): The 10/im emission profile along the ,v-coordinate at y = 0, 
normalized to its central value. (b3): Same as (b2), but vertically along the y-coordinate at x = 0. 



it is flat with ^ = 0. This leads to a strongly constrained torus 
thickness Y = 20, in sharp contrast with SST1721-I-6012 where 
Y is the least-well constrained parameter. All 28 models in 
the selected sample have the same value of Y, although this 
probably reflects our discrete sampling of parameter space; 
there could be a small range around Y = 20, but 10 and 30 are 
clearly excluded. The cloud optical depth is well-constrained 
at Tv « 20 - 30. On the other hand, while well constrained 
for SST1721-I-6012, A^o is the least well constrained parameter 
here, with a likely range of 2-9 clouds. 

As noted above, flat-top emission features arise from ab- 
sorption by a single cloud with Ty ^ 20. As is evident 
from Table 3, all accepted models have Ty ^ 20-30, but A^o 
is largely unconstrained. However, the average number of 
clouds along the line of sight to the AGN (see Equation 3) 
falls within the narrow range < A^(/) < 2.5 for all accepted 
models. To a certain degree, A^(/) is a good proxy for the 
typical number of clouds along lines of sight that pass close 
to the dust sublimation radius, where the 10/im emission is 
originating. In panel (al) of Figure 7 we show the number of 
clouds along all lines of sight through the best-fit torus model. 
Panel (a2) shows a one-dimensional cut through the image, 
providing the profile of the number of clouds per ray along 
the x-coordinate at y = 0. Due to the axial symmetry of the 
torus, this profile is symmetric with respect to x = 0, irrespec- 
tive of the viewing angle. The signature of the central cav- 
ity is clearly visible in this profile: the cloud column reaches 
a minimum of 2.4 at the center and stays close to this level 
for all x/Rj < 1. It reaches a maximum of w 6 clouds along 
rays roughly lORd away from the AGN. Panel (a3) shows the 
corresponding profile in the vertical direction at x = 0. The 
symmetry of this profile around y = Q again reflects the axial 
symmetry, which ensures equal path lengths through the torus 
above and below the central line of sight. 



While the number of clouds along two lines of sight dis- 
placed symmetrically from the center is equal, the illumina- 
tion patterns of individual clouds as seen by the observer can 
differ for the two, depending on the position angle in the plane 
of the sky. In panel (bl) we plot the two-dimensional dis- 
tribution of the model 10/im emission for the central region 
with size lORd x lORj. Roughly 50% of the flux is detected 
within the inner 5Rd radius, and 70% of that fraction comes 
from the image upper half. This radiation originates from re- 
gions on the far inner face of the torus; no emission originates 
from the near side, where the observer faces the dark sides of 
the clouds. Similar to panels (a2) and (a3), we plot in pan- 
els (b2) and (b3) profiles of the 10/im emission along the x 
and y directions. As expected, the horizontal profile in panel 
(b2) is symmetrical, clearly displaying the dust-free cavity at 
its center On the other hand, the shape of the vertical pro- 
file in (b3) reveals the asymmetry between the emission in the 
upper and lower halves. Despite equal cloud columns along 
viewing lines above and below the image center, the emission 
is not equal, owing to the strong anisotropy of single cloud 
emission. The 10/tm emission originates from hot, bright sur- 
faces of clouds located on the torus far inner face. Clouds in 
the torus near side, which show their dark, cooler faces, only 
absorb the 10/im photons that were emitted on the torus far 
side. 

Most clumpy torus models do not produce the apparent shift 
in peak emission. The shifts occur predominantly in models 
that have a small Ty (< 20). Significantly, Ty ^ 20 mod- 
els are also the ones producing the most prominent 10/im 
emission features across the likely range of Ty. As is evident 
from Fig. 16 in N08b, the emission feature strength decreases 
monotonically as Ty increases up to Ty ~70; in some cases 
the feature even switches to absorption for pole-on viewing. 
Therefore low-Ty models stand out in their feature strength 
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and it is reasonable that such sources would be preferentially 
selected in observations that looked to identify the lOfim sili- 
cate emission feature in AGN. Finally, it should be noted that 
the absorption coefficients widely used in the literature do not 
have their peaks at 9.8 /im. In the tabulation of Draine (2003), 
the feature peaks at 9.48 fim instead. The "cold" silicate dust 
of OHM, which is the one used here, has its peak at 10.0 /im. 
The effect on the present discussion is insignificant. 

4. SPECTRAL PROPERTIES OF CLUMPY MODELS 

In addition to the detailed fitting of the SST1721+6012 and 
PG121 1+143 data, we investigated some other properties of 
the 10/im feature, comparing observations with general prop- 
erties displayed by the model database. 

Hao et al. (2007) present a large compilation of Spitzer mid- 
IR spectra. Although a loosely defined sample, it is the largest 
gathered thus far For each source they measure the 10/im 
silicate feature strength ^lo from 



>io ■ 



In 



^cont(Alo) 



(4) 



where F is the measured flux, Fcom is a continuum constructed 
underneath the 10/im and 18/im silicate features (see Sirocky 
et al. 2008, for details; see also §3), and Aio is the peak wave- 
length of the feature strength; emission features have a posi- 
tive ^lo, absorption features a negative one.^ It may be noted 
that the specific prescription of continuum construction mod- 
ifies, and can even reverse, the relative strengths of the 10/im 
and 18/im features, as is also apparent from Figure 6; the ratio 
of the two strengths is an important indicator of dust optical 
properties (Sirocky et al. 2008). 

Removing all ULIRGs, the Hao et al. (2007) sample con- 
tains 21 QSOs, 38 Seyfeit 1 and 39 Seyfert 2 galaxies. The 
top panel of Figure 8 shows the histograms of the feature 
strengths for the three groups. The figure other panels show 
results from recent studies, which produced additional compi- 
lations of feature strengths: Thompson et al. (2009) compared 
a sample of Seyfert 1 galaxies with quasars, Wu et al. (2009) 
and Gallimore et al. (in preparation) analyzed Seyfert galax- 
ies, both type 1 and 2, from the 12/im Galaxy Sample (Rush 
et al. 1993). While Wu et al. adopt the original source clas- 
sification of Rush et al., Gallimore et al. establish a different 
source type in several cases, based on work published else- 
where. We employ the latter classification in both panels (c) 
and (d), dispensing with all sources re-classified as LINERs or 
HII (star-forming) galaxies and confirming the Wu et al. sug- 
gestion that the re-classification of several sources has little 
effect on the statistical results of 5io measurements. 

In addition to the torus emission, the infrared radiation of 
many active galaxies contains a starburst contribution whose 
fractional strength varies from source to source (e.g., Netzer 
et al. 2007). Removing the starburst component by subtract- 
ing a suitable template and leaving no PAH residuals is thus 
an important preliminary step in the detailed SED analysis 
of many individual AGN (see, e.g., Mor, Netzer, & Elitzur 
2009). Note, however, that the two sources analyzed here in 
detail show no signs of ongoing star formation, either in the 
form of PAH emission or far-IR (FIR) emission. While PAH 



' The feature strengths of SST1721+6012 are Sio ■ 
(for the feature). Both are uncertain to within ' 



= 0.26 and Sig = 0.34 
' it 0. 1. For the model 

shown in Figure 6, the feature strengths determined from the blue curve are 
5io = 0.26 and 5i8 = 0.21 (at 18.0 ^im), and Sio = 0.33 and Sig = 0.21 (at 
17.5 /xm) from the red curve. 
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Fig. 8. — Distributions of the lO^m silicate feature strength S\q (Equa- 
tion 4) in several AGN samples. The bin size is 0.1. All ULIRGs present in 
the original samples have been removed. Measurements for QSOs are shown 
with gray bars, Seyfert 1 s with dashed lines, and Seyfert 2s with solid lines. 
The number of sources of different type is given in parentheses in the leg- 
end, (a): Spitzer sample by Hao et al. (2007). (b): Archival sample of type 1 
sources by Thompson et al. (2009). Note the different scale, (c): Seyfert 
sources from the 12/im Galaxy Sample, presented by Wu et al. (2009), and 
(d) re-analyzed by Galhmore et al. (in preparation). Panel (c) contains only 
sources also present in (d), and the source classification in both panels is 
adopted from the latter. 



emission can contaminate the 10/im region in some individ- 
ual spectra, its overall impact on the averages of large sam- 
ples seems minimal. Netzer et al. (2007) subtract a starburst 
template from the average spectra of AGN with and without 
strong FIR detections and find that the MIR regions are hardly 
affected by this subtraction in either case. In particular, their 
Figure 6 shows that the strength of the 10/im silicate feature 
barely changes. The analyses by Wu et al. (2009) and Gal- 
limore et al. (in preparation) of the same data set provide an 
even stronger evidence: The former ignores the potential star- 
burst contribution while the latter includes a PAH component, 
handled with the Pahfit tool (Smith et al. 2007). In spite 
of this difference, the histograms in panels (c) and (d) of Fig- 
ure 8 are quite similar, showing comparable lower and slightly 
increased upper limits on ^lo and an overall shape that is es- 
sentially the same. 

Comparison of the histograms for type 1 sources in the pan- 
els of Figure 8 shows that in moving from Seyfert to quasar 
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luminosities the 10/im feature shifts to enhanced emission. 
This trend was noted earlier in Nenkova et al. (2008b; see 
§6.4), and the analysis here verifies this suggestion, giving it 
quantitative evidence. Nenkova et al. point out that the most 
likely explanation is that the number of clouds along radial 
rays is smaller in quasars than in Seyferts. 

Grouping together the QSOs and Seyfert Is of the Hao 
et al. (2007) sample, the top panel of Figure 9 shows the his- 
tograms and Table 4 lists the statistical indicators of the Sio 
distributions in type 1 and 2 sources. Most sources exhibit 
rather small absolute values of ^iq. The histogram of type 1 
sources is clearly shifted toward emission in comparison with 
type 2. Although the Hao et al. (2007) sources do not consti- 
tute a complete sample, the selection criteria were unrelated 
to the silicate feature. The derived histograms can thus be rea- 
sonably considered representative of the differences between 
types 1 and 2. 

Our clumpy torus models should produce similar his- 
tograms if they bear a resemblance to the IR emission from 
AGN. Such a comparison presents two fundamental difficul- 
ties. First, the assignment of a given clumpy model to type 1 
or 2 is not deterministic — only a probability can be assigned. 
We handle this problem by dividing the models according to 
the probability Pesc for an unobscured view of the AGN. The 
collection of models with Pesc > 0.5 can be expected to resem- 
ble the behavior of the type 1 population, those with Pesc < 0.5 
type 2. The second problem is that the actual distribution of 
parameter values is unknown. Since we do not have any han- 
dle on these distributions, we decided to test the adequacy 
of histograms produced by a uniform sampling of the model 
parameters within the bounds deduced in N08b: < q < 3, 
No < 15, 30 < TV < 100, 15° < o- < 60°, and 10 < F < 100. 
Since q was sampled here more thoroughly than in NOSb, we 
use the full range listed in Table 1. The parameters were 
sampled in steps of 0.5, 1, 10, 5, 10, 10 for q, No, Ty, a, 
y, /, respectively. The bottom panel of Figure 9 shows the 
histograms of S\o for all database models selected by these 
criteria. These distributions resemble those of the observa- 
tional sample, as is also evident from their statistical proper- 
ties listed in Table 4. The lo-, and 2<Js ranges of ^lo given in 
the table contain sources and models within 1 and 2 standard 
deviations CTj from the mean of each of the two distributions. 
Leading to exclusion of only few sources and models at the 
very ends of the distributions, this additional selection has the 
effect of a much more meaningful agreement of the two dis- 
tribution widths, not spoiled by rare outliers. At the level, 
rejected sources are just 3 type 1 and 2 Seyfert 2s, and among 
the sample of models only 4% of those with Pesc > 0.5 and 6% 
with Pesc < 0.5 are excluded due to this criterion. 

Although the choice of uniform sampling of the model 
database is arbitrary, it produces reasonable results. The rea- 
son is that, as noted already in N08a and NOSb, clumpy mod- 
els never produce very deep absorption features, in agree- 
ment with observations. This limited range is reflected in the 
histograms for any reasonable criteria used for model selec- 
tion from the database. The other main characteristic of the 
observed histograms is the separation between type 1 and 2 
sources, and this, too, is reproduced reasonably well by the 
models. 

5. SUMMARY AND DISCUSSION 

Spitzer IR observations of AGNs have increased signifi- 
cantly the number and quality of SEDs for these objects and 
produced some puzzling results, especially with regard to the 
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Fig. 9. — Distributions of the silicate feature strength Sio- The bin 

size is 0.1 and each histogram is normalized to unit area. Top: Data from 
the Hao et al. (2007) sample. Dotted line shows type 1 sources (QSO and 
Seyfert 1 combined), solid line Seyfert 2s. Bottom: Histograms for 840,000 
Clumpy models whose parameters most likely correspond to physical values 
(see text). Models with escape probability P^sc > 0.5 (likely type 1 source in 
a clumpy torus) are shown as dotted line, those with Pcsc < 0.5 (likely type 2) 
as solid line. For statistical properties of all samples see Table 4. 



TABLE 4 
Sio STATISTICS 





Hao et al. 


(2007) 


Clumpy 


source type 


QSO + Syl 


Sy2 


Pesc > 0.5 


fesc < 0.5 


sample size 


59 


39 


340,000 


500,000 


mean 


0.03 


-0.46 


0.15 


-0.33 


median 


0.12 


-0.34 


0.12 


-0.32 




0.36 


0.40 


0.29 


0.38 


Ictj -range" 


-0.30, 0.37 


-0.83, -0.06 


-0.14, 0.44 


-0.71,0.04 


2(T., -range'' 


-0.44, 0.40 


-1.20, 0.29 


-0.42, 0.73 


-1.09, 0.42 



Note. — Samples of AGN and of Clumpy models as in Figure 9. 
" Ranges of 1 standard deviation CTj from a sample's mean value. 
The 2-standard deviations range. 



10/im silicate feature. These include (1) detection of the fea- 
ture in emission in type 2 sources, (2) emission features with 
broad, flat-topped peaks shifted toward long wavelengths in 
several type 1 sources, and (3) absence of any deeply absorbed 
features. None of these observations can be satisfactorily ex- 
plained with smooth density torus models. 

Here we have shown that clumpy torus models provide rea- 
sonable explanations for all three puzzles. To that end we have 
fitted the Spitzer SEDs of two very different sources with our 
Clumpy models. One source, SST1721H-6012, is the first 
type 2 QSO to show a clear 10/im emission feature. Our anal- 
ysis provides a reasonable fit of the SED with a model that 
shows the feature in emission. In contrast with smooth den- 
sity models, where the AGN is either obscured or visible, our 
model produces a small obscuration probability, Pobsc = 27%, 
for this type 2 source. This relatively low probability may 
explain why SST1721H-6012 is the only source among more 
than 20 type 2 QSOs with measured SEDs (see, e.g., Polletta 
et al. 2008) to show a clear 10/im emission feature. 

Addressing the second puzzle, PG1211H-143 is one of the 
first QSOs to display the 10/zm silicate feature in emission. 
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a feature that is unexpectedly broad and apparently shifted to 
longer wavelengths. The original attempts to explain these 
properties invoked non-standard chemical dust composition. 
Our modeling shows that the shifts are only apparent and re- 
sult from the flattening of the feature peak by radiative transfer 
in clumpy media. The feature is well reproduced by clumpy 
models with standard dust. The third observational puzzle, 
lack of deep 10/im absorption features in any AGN, has al- 
ready been shown to be a signature of clumpy dust distri- 
butions (Nenkova et al. 2002; Levenson et al. 2007; Sirocky 
et al. 2008; Nenkova et al. 2008 a,b). Here we go a step fur- 
ther and produce the histogram of 10/im feature strength for a 
large sample of AGN CLUMPY models. The result is in good 
qualitative agreement with the sample observed by Hao et al. 
(2007). In particular, the median values of both type 1 and 
type 2 observed distributions and their widths are well repro- 
duced by the model database. 

The IR SED generally does not constrain very tightly the 
properties of dusty sources — the large degeneracy of the ra- 
diative transfer problem for heated dust is well known (e.g., 
Vinkovic et al. 2003). In the present case, the problem is fur- 
ther exacerbated by the clumpy nature of the dust distribution 
and the non-spherical geometry. Our model database contains 
close to 5 million entries, and although the fitting procedure 
eliminates most of them, many produce reasonable agreement 
with the observations. In the case of SST1721H-6012, close to 
1,700 models with very different parameters deviate by no 
more than 10% from the best-fit model. And while the much 
higher quality of data in PG1211H-143 greatly reduces the 
number of acceptable models, there are still 28 different ones 
that are practically indistinguishable in the quality of their fits. 
In the face of this degeneracy, we have developed a statistical 
approach to assess the meaningfulness of the various torus pa- 
rameters derived from the fits. We find that some parameters 
are well constrained in each case, while others are not. In 



both sources the power law of the radial distribution (q) and 
the optical depth of a single cloud (tv) are well constrained, 
while the torus viewing angle (/) and its angular thickness (cr) 
are not. Both the cloud number (A^o) and radial thickness (Y) 
are well constrained in only one of the sources, a different one 
in each case. Asensio Ramos & Ramos Almeida (2009) have 
recently developed a different, novel approach to tackle the 
degeneracy problem. They interpolate the CLUMPY SEDs by 
means of an artificial neural network function, allowing them 
to study the parameter distributions as if they were contin- 
uous, and employ Bayesian inference to determine the most 
likely set of parameters. Applying this method to a selection 
of sources, Ramos Almeida et al. (2009) find that the princi- 
pal ability to constrain different CLUMPY parameters strongly 
depends on the individual source. We have already begun an 
extensive comparison of the two approaches and will report 
our findings elsewhere. 

Although the SED alone is generally insufficient for de- 
termining all the torus parameters with certainty, the success 
in resolving outstanding puzzling behavior of the 10/im fea- 
ture in AGN is encouraging and enhances confidence in the 
clumpy torus paradigm. 
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